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High-velocity galactic outflows, driven by intense bursts of star formation and black hole accretion, 
are invoked by current theories of galaxy formation 1 to terminate star formation in the most massive 
galaxies and to deposit heavy elements in the intergalactic medium. From existing observational 
evidence 2 3 on high-redshift galaxies, it is unclear whether such outflows are localized to regions of 
intense star formation just a few kiloparsecs in extent, or whether they instead have a significant 
impact on the entire galaxy and its surroundings. Here we present two-dimensional spectroscopy of 
a star-forming galaxy 4 at redshift z=3.09 (seen 11.5 gigayears ago, when the Universe was 20 per 
cent of its current age): its spatially extended Lya line emission appears to be absorbed by HI in a 
foreground screen covering the entire galaxy, with a lateral extent of at least lOOkpc and remarkable 
velocity coherence. It was plausibly ejected from the galaxy during a starburst several 10 8 years 
earlier and has subsequently swept up gas from the surrounding intergalactic medium and cooled. 
This demonstrates the galaxy-wide impact of high-redshift superwinds. 



The formation of galaxies requires gas to cool in haloes 
of dark matter that collapse under gravity from the expan- 
sion of the Universe. However, cooling alone overproduces 
bright galaxies at the present day, so models incorporate 
thermal conduction, photoionization and galaxy merg- 
ing, together with additional 'feedback' 1 in the form of 
galactic-scale outflows. The latter are powered by super- 
novae and massive stellar winds, or by relativistic winds 
and jets resulting from gas accretion onto supermassive 
black holes. Although starburst superwinds have been 
studied 5 in local dwarf galaxies (such as M82), observa- 
tional evidence for their counterparts in young massive 
galaxies at high redshift has been less direct. 

In Lyman-break galaxies (LBGs) - high-redshift galax- 
ies with moderate masses and star-formation rates - spec- 
troscopic evidence points to powerful outflows 2 ' 3 of inter- 
stellar gas, with absorption lines being blueshifted by sev- 
eral hundred kilometres per second from the galaxy's sys- 
temic velocity. The gravitationally lensed LBG MS 15 12- 
cB58 at redshift z=2.73 is a clear example 3 : absorbing 
gas is outflowing at ~ 255 km s -1 at a rate exceeding the 
star-formation rate, and although it covers the entire star- 
forming region, this shows only that the radius of any shell 
exceeds ~ 1 kpc; it is not known whether the outflows 
are collimated and localized to star- forming regions, or 
whether they are galaxy-wide. The latter is suggested by 
the observed decrease in the intergalactic medium (IGM) 
HI opacity in background quasar sightlines that are close 



in projection to LBGs 6 (within 0.7 Mpc co-moving). This 
may indicate that LBGs drive superwinds at ~ 600 km s _1 
for several 10 s yr, or it might be a statistical fluctuation 
since there are few LBG-quasar pairs at the smallest sep- 
arations. 

To characterize such outflows better via absorption 
studies, a background light source is needed with a spa- 
tial extent somewhat larger than an LBG stellar body 
or a quasar sightline. Such a source is provided by the 
recently discovered Lya-emitting blobs (hereafter LABs), 
associated with LBGs in the SSA22 protocluster at red- 
shift z=3.09, a structure which is likely to evolve into 
a rich cluster of galaxies. With sizes ~ 100 kpc and 
Lya luminosities ~ 10 44 erg s _1 , the LABs resemble the 
emission-line haloes of high-redshift radio galaxies, but 
their radio emission is much fainter than in typical high- 
redshift radio galaxies 7 . Models for LAB Lya emission in- 
clude superwind ejecta 8 and cooling radiation from grav- 
itationally heated gas 9 . Both LABs are submillimetre 
sources 10 , with LAB-1 the more powerful and consistent 
with ~ 1000 M Q yr -1 of star formation. LAB-2 has an 
X-ray detection, suggestive of an obscured active galactic 
nucleus 11 , whereas LAB-1 does not. 

We recently observed the LAB Lya haloes using 
integral-field spectroscopy, which, unlike conventional slit 
spectroscopy, gathers spatially resolved spectra over a 
two-dimensional area. Such information is essential for 
a complete picture of morphologically complex objects 
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Figure 1. Part of the data cube obtained from integral field 
spectroscopy of the z=3.09 galaxy LAB-2, showing the Lya 
emission line. The absorption feature, interpreted as due to 
foreground HI, in indicated. A data cube is a three-dimensional 
structure for spectroscopic data from integral field spectro- 
graphs. Here the long axis corresponds to the dispersion axis 
and the other two axes represent spatial dimensions on the 
sky, that is, the data cube stores a spectrum for each point in 
a two-dimensional region of sky. With visualization software, 
the data-cube representation can provide useful qualitative in- 
sights into complex data sets. A, wavelength. 

like the LABs. We used the SAURON Integral Field 
Spectrograph 12 on the William Herschel Telescope, pro- 
viding moderate resolution spectroscopy (~ 4A full-width 
at half-maximum, FWHM) over a 41 x 31 arcsec area, 
sampled with 0.95-arcsec lenslets. LAB-2 was observed 
for 15 h (in thirty 30-min sub-exposures), but unlike LAB- 
l 13 exhibits simple structure in Lya. At most positions 
the emission line is double-peaked, as shown in Fig. 1 
and in the line profiles of individual regions (Figs. 2 
and 3). The position of the profiles central minimum 
varies across the galaxy by less than 60 km s _1 in its 
rest-frame. The observing strategy of dithering between 
sub-exposures ensures that this feature is not a charge- 
coupled device (CCD) detector defect, nor is it due to a 
night-sky emission feature. The simplest interpretation 
is that part of the Lya emission is scattered out of our 
line of sight by neutral hydrogen in a coherent foreground 
screen. 

The line profiles were each modelled with a gaussian 
emission line and a superimposed Voigt profile absorber 
(all convolved with the instrumental response), as in Fig. 
3. The HI column densities are Nhi ~ 10 19 cm~ 2 , and the 
relative amplitudes of the red and blue line wings reflect 
velocity differences between the emission and absorption 
components. The absorber is typically blueshifted from 
the underlying peak by less than 250 km s _1 ; the emis- 
sion shows a range in absolute velocity of ~ 290 km s _1 , 
with individual linewidths of ~ 1000 km s^ 1 FWHM. At 
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Figure 2. A reconstructed image of the z=3.09 galaxy LAB-2 
in the Lya emission line. This shows the extent of its gaseous 
halo (10 arcsec on the sky equates to a physical distance 
of 76 kpc in this galaxy with the cosmological parameters 
H = 70 km s" 1 Mpc~\ M = 0.3 and ft A = 0.7). The im- 
age was derived from observations with the SAURON integral 
field spectrograph, by collapsing the data cube shown in Fig. 
1 along the wavelength direction. The labels indicate regions 
for which one-dimensional Lya emission line profiles have been 
extracted, as shown in Fig. 3. The underlying Lya emission is 
likely to be from gas in the halo of the associated Lyman-break 
galaxy 4 (at position C), arising from a combination of super- 
wind ejecta 8 , cooling radiation 9 , and gas photoionized by the 
obscured active nucleus 11 . 

position B, no absorber is detected, perhaps because the 
underlying emission is more blueshifted than elsewhere. 

An absorbing shell is predicted to form when a 
starburst-heated hot gas bubble becomes over-pressurized 
relative to the interstellar medium and hence breaks out 
of the galactic disk, accelerates, and fragments as a re- 
sult of Rayleigh- Taylor instabilities. Hot gas then es- 
capes into the halo and forms a second shell of swept-up 
IGM. Evolutionary models for the Lya emission from such 
superwinds - which are consistent with Hubble Space 
Telescope observations of local starburst galaxies 15 - sug- 
gest that the LAB-2 absorber is in a late phase, where 
the shell has cooled and slowed sufficiently to absorb the 
underlying emission. Figure 4 illustrates the model. 

Because the absorbing shell's transverse extent is at 
least 100 kpc, it is probably a comparable distance in front 
of the Lya source (although the constant absorber velocity 
implies that, if spherical, the radius must be considerably 
larger). At a velocity of 1000 km s _1 , it takes 10 8 yr to 
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Figure 3. Lya emission line profiles across the galaxy LAB-2 
at z—3.09. A profile is shown for each region labelled in Fig. 
2. Thick solid lines are fits to a model where the intrinsic Lya 
emission (dashed lines) is partially absorbed by foreground HI. 
The measured HI column densities are ~ 10 19 cm~ 2 , although 
the modest spectral resolution implies la uncertainties of at 
least a factor of ~ 10, so the column density may not be 
uniform across the galaxy. The widths of the absorbers (ex- 
pressed as the gaussian b-parameter in the Voigt profile) are 
~ 20 km s~ , which for thermal broadening implies a kinetic 
temperature ~ 10 4 K (as expected on photoionization grounds 
and consistent with the lack of bulk motion across the shell). 
Note, however, that similar Lya profiles can result from a Lya 
source in an HI slab of much higher optical depth 19 ; the twin 
peaks then arise because the escaping photons scatter into the 
line wings where the optical depth is lower. This would not, 
however, explain why the apparent absorption feature is con- 
stant in wavelength while the fitted emission velocity varies by 
almost 300 km s _1 . 

travel 100 kpc. Although 1000 km s -1 is the expected 
level of the initial velocity, it has since slowed owing to 
mass loading from the IGM and work against the galaxy's 
gravity; this implies an age of several 10 8 yr. 

We now consider the cooling of the shocked gas. Gas 
behind a shock of speed v s h oc k is heated to T = 1.4 x 10 
(vshock/1000 km s~ x ) 2 K, from which the radiative cool- 
ing time is i coo i = l-5fcT/4nh a ioA, where A is the cool- 
ing function 16 (> 5 x 10~ 23 erg cm 3 s _1 for 10 6 ~ 7 K 
and solar abundances) and nhaio is the surrounding gas 
density (the factor 4 reflects the density increase behind 
a strong shock). Together, these imply that i coo i < 
4.6 x 10 5 yr (v shock /1000 km s- 1 ) 2 /(n ha i / Cm - 3 ). Moving 
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Figure 4. An illustration of our model for the Lya emis- 
sion and absorption in the redshift z=3.09 galaxy LAB-2. The 
Lya emission is produced by gas with a velocity dispersion 
~ 1000 km s _1 and bulk motions ~ 290 km s _1 . This emis- 
sion is then absorbed/scattered by HI in a foreground shell. 
The shell's velocity coherence across the source suggests that 
it is over 100 kpc from galaxy, having been swept up from 
the surrounding medium by a starburst-driven outflow sev- 
eral 10 8 yr earlier. The possibility that the absorption is 
not directly associated with LAB-2, but part of a larger fore- 
ground absorber which just happens to cover it, is unlikely 
since this would imply an overall size much larger than ab- 
sorbers of comparable Nhi in quasar sightlines 20 : the LAB-2 
absorber column density is between that of a Lyman-limit 
system (diameter ~ 80 kpc) and that of a damped Lya ab- 
sorber (~ 30 kpc). Such sizes are consistent with recent 
modelling 21 and cosmological simulations 22 . Moreover, it has 
been shown 23 that all damped Lya absorbers at redshift z ~ 3 
with Nhi > 2 x 10 20 cm~ 2 can be explained by LBG galactic 
winds, if they are driven out to ~ 27 kpc. 



out from the galaxy into the IGM at 100 kpc, the den- 
sity nhaio will range from a few 10~ 3 to a few 10 -4 cm~ 3 , 
so for shock velocities of several hundred kilometres per 
second, the shocked gas can cool within the expansion 
timescale of a few 10 8 yr. The gas is expected 17 to cool 
to ~ 10 4 K and come into photoionization equilibrium 
with the metagalactic ultraviolet radiation field, with an 
HI fraction ~ 0.1 — 0.001. Indeed, for the above n halo 
values, the total column density swept up out to 100 kpc 
could be up to 10 21 cm~ 2 , compared with the measured 
Nhi — 10 19 cm~ 2 . The implied mass of the shell, if it 
covers an area 100 kpc square, could be almost 10 11 M©. 
For a velocity of 250 km s _1 the kinetic energy of such 
a shell is 6.3 x 10 58 erg, which we can compare with the 
energetics of supernovae: the LBG associated with LAB- 
2 (M14) has a near-infrared Ks-band magnitude 4 of 21.7, 
which implies (from galaxy evolution models 18 ) a baryonic 
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mass of ~ 10 M Q , mostly in stars. For a normal stellar 
initial mass function, supernovae provide ~ 10 49 erg per 
solar mass of stars 1 , so ~ 10 60 erg is available and the 
superwind model appears energetically feasible. 

We thus have a consistent picture in which the Lya 
absorber is a highly ionized shell of gas swept up from the 
IGM over several 10 8 yr since the end of the starburst, 
having slowed to within a few hundred kilometres per sec- 
ond of the galaxy's systemic velocity. This demonstrates 
that superwinds are a galaxy-wide phenomenon, whose 
impact is not confined to small-scale regions of current 
star formation, as previous observations might have sug- 
gested. This underscores their importance as a primary 
feedback agent in galaxy formation. 
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